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Abstract. We present results on dynamical processes that exhibit a stretched 
exponential relaxation. When the relaxation is a result of two competing expo- 
nential processes, the size of the system, although macroscopic, play a dominant 
role. There exist a crossover time t x that depends logarithmically on the size of 
the system, above which, the relaxation changes from a stretched exponential to a 
simple exponential decay. The decay rate also depends logarithmically on the size 
of the system. The results are relevant to large-scale Monte-Carlo simulations and 
should be amenable to experiments in low-dimensional macroscopic systems and 
mesoscopic systems. 



Many relaxational processes in macroscopic systems are characterized by a 
relaxation function Q(t) that exhibits a stretched exponential behavior, 

Q(i)~Q(0)exp[-(t/r) /3 ], (1) 



where < f3 < 1. Examples include viscoelastic relaxation 
ation (2), glassy relaxations f§,4j|, relaxation in polymers 
decay in trapping processes [S|. Many more examples j£ 
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suggest that ( jj is common to a very wide range of phenomena and macro- 
scopic materials. 

The origin of the stretched exponential is not always clear. In many cases 
it is assumed to be the result of a competition between two exponential pro- 
cesses. In some cases, e.g., trapping processes at long times, this assumption 
is well established, while in others, such as relaxation in glassy materials, this 
assumption has been controversially discussed |l(||l7| and alternative models 
have been also suggested 0,^8 19 2(J 



We have recently investigated the occurence of stretched exponential be- 
havior in finite systems, in cases where the relaxation arises due to two com- 
peting exponential processes [ pT| 

We have found that: (a) the size of the system, although macroscopic, 
plays a dominant role in the relaxation time pattern, leading to an expo- 
nential decay sufficiently at long times; (b) the crossover time, t x , to the 
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exponential depends logarithmically on the system size; (c) the rate of the 
exponential decay also depends logarithmically on the system size, and (d) 
in the special examples of the trapping and the hierarchically constrained 
dynamics models the exponential relaxation may enter before the stretched 
exponential is reached. These results are of relevance to experiments in con- 
fined systems, mesoscopic systems and to Monte-Carlo simulations. Our the- 
oretical predictions on the finite size effects can serve as an experimental test 
for identifying the origin of the mechanism leading to stretched exponential 
decay. 

We assume that the relaxation function of the whole system can be rep- 
resented by an integration over all possible states n, namely, 



Here, #(n) is the probability that state n is occupied and Q(n,t) is the 
dynamic relaxation of the n-th state. 

Usually, in the case of a stretched exponential behavior, <£(n) is assumed 
to behave as <P(n) ~ exp(— an a ), while Q(n,t) decays exponentially with 
time as Q(n,t) ~ exp(— bt/n 1 ). A number of dynamical models that yield a 
stretched exponential decay can be formulated in terms of Eq.(^). These in- 
clude the long-time behavior in the trapping problem Q , the target problem 
[ pp| , direct energy transfer ]2p| , trapping of nonidcntical interacting parti- 
cles p3fl ) hierarchically constrained dynamics [jl6|, models for relaxation in 
microenulsions and molecular glasses and others. 

We now concentrate on three examples: (i) A particle diffusing in a d- 
dimensional system with randomly distributed static traps, where we are 
interested in the survival probability Q(t) of a particle. Here the state n rep- 
resents a particle in a trap-free region of linear size n; <P(n) is the probability 
for the occurance of a size n trap-free region, and Q{n, t) is the survival prob- 
ability of the particle in this region II] . The exponent a is the dimension d 
of the system, and 7 = 2 due to the diffusional motion, (ii) A linear system 
(chain) along which two types of particles {A and B) are diffusing and in- 
teracting via hard core interaction. However, only type A can be trapped by 
static traps which are randomly distributed along the chain. Here, Q(t) is the 
survival probability of particles of type A, <P(n) is the probability that a free 
trap region of size n occurs, and Q(n,t) is the survival probability of a type 
A particle to survive in this region. The exponent a is the dimension of the 
system a — 1 and 7 = 4 is due to diffusion in the presence of hardcore iter- 
ations p5[ . (iii) Hierarchically constrained dynamics, a model that has been 
proposed to account for glassy relaxation jl6| . This model assumes that the 
relaxation of level n populated by spins, occurs in stages, and the constraint 
imposed by a faster degree of freedom must relax before a slower degree of 
freedom can relax. This implies that the time scale of relaxation in one level 
is subordinated to the relaxation below. A possible realization considered in 
[ fl6| and here is a system with a discrete series of levels where the relaxation 
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time of level n is r n ~ n 7 (corresponding to the exponential form of Q(n,t) 
in (|2]) ) , and the weight factor of level n, is <£(n) ~ e~ an JL2| , corresponding 
to a = 1. The first exponential in (|^) is accordingly the probability to occupy 
level n and the second exponential represents the decay of that level. 

The evaluation of the long time behavior of the integral in (||) is performed 
using the method of steepest descent. The main contribution to the integral 
arises from the maximum of the integrand in (Q) , which is obtained from the 
minimum of the function, —an a — bt/rf 1 , appearing in the exponent. This 
yields that the main contribution to (2) comes from 

n* = ( 7 6i/aa) 1/(a+7) , (3) 

leading to [l] with (3 = a /(a + 7) < 1, and r = (a/b^a-t^i-y/ij + a)) 1+ ~i/ a . 

However, as shown below, these arguments are valid only in the thermo- 
dynamic limit where the system size is infinite. For a finite number N of 
traps (in the trapping system) or a finite system with a finite number N of 
spins (in the hierarchical constraint system) the relaxation function depends 
explicitly on N . Since our discussion is quite general for systems described by 
(^|), in what follows we refer below to traps and spins in the above examples 
as elements. 

For a single finite system consisting of N elements, the relaxation function 
Q{t) represents an average quantity over the N elements, 

Q(t) = ^^m(n)Q(M), (4) 
{«} 

where the sum is over all possible states n and m(n) is the number of elements 
at state n, with m ( n ) — N. Since the sum in (^) is over exponential 

functions, the value of Q(t) will fluctuate for different sets of N. There will 
be a distribution of Q(t), and we are interested in the typical Q(t), which is 
around the peak of this distribution. 

In the thermodynamic limit N — > 00, all states n are occupied, m(n)/N 
can be identified with <P(n) and (Q) follows. For N finite, in contrast, there 
exists a characteristic "maximum" state n = n max (./V), and this n max should 
replace the upper limit (00) in (^|), 

Q(t) = I ^ $(n)Q{n,t)dn. (5) 
Jo 

To estimate how n max depends on N, we note that the typical number of 
states n in a sample of N elements is Z(n) = N<P(n) = 7Vexp(— an a ). States 
with Z(n) <C 1 will not occur in a typical system of N elements, and this 
yields 

hi^V\ 1/Q 

(6) 
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If n* <C rimax, the upper limit in can be approximated by infinity and 
thus leads to (^). However, if n* 3> n max the main contribution to (^) will 
not be from the maximum of the integrand, which is outside the range of 
integration, but from n max . Thus, for n* 3> n max we expect 

Q(t) S Q(0)e- b */"-* (7) 

where the time constant of the relaxation, nX ax , scales as (lniV) 7 / Q . The 
crossover time from a stretched exponential (|l[) to an exponential (fjj) can be 
estimated from the condition n* — ro max , from which follows 

aa flnN\ 1+l/a 

*x = -T — ■ (8) 



7& \ a 

The striking point in (^) is the logarithmic dependence on N, which puts 
t x in the range of observable time scales measurable in mesoscopic and even 
macroscopic systems. Indeed, the corresponding relaxation value Q(t x ) scales 
as 

Q{t y )~N- a '\ (9) 

independent of the microscopic parameters a and b. For the above three cases 
we find: (i) In the case of the trapping relaxation mechanism where a = d 
and 7 = 2 we obtain, 

Q(t x )/Q(0)~N- d / 2 . (10) 

(ii) In the non identical case a = 1 and 7 = 1 and thus 

Q(U/Q(o) ~ A^ 1/4 . (11) 

(iii) In the hierarchical constraint dynamics 

O(*x)/O(0) -N- 1 ^. (12) 

It is known [8e,23] that in both examples, for an infinite system, the stretched 
exponential behavior of (Q) sets in only at very long times.. Thus we expect 
that in the finite system, the crossover will mask the stretched-exponential 
pattern. 

To test our analytical approach, we performed new Monte Carlo sim- 
ulations on two cases (i) and (iii), the trapping model (case (i)) and the 
hierarchical constraint model (case (iii)). In the trapping model, we consider 
one and two dimensional systems with a fixed concentration c = 0.5 of ran- 
domly distributed traps, and vary the size N/c of the system. We calculated 
numerically the survival probability Q(t) of a particle as a function of t and 
N. In the hierarchical model we have chosen t„ ~ n i.e., 7=1. We calculated 
the relaxation function for system sizes varying from N — 10 2 to N = 10 5 . 

As mentioned earlier, the relaxation function fluctuates for different sets 
of N. For obtaining the typical behavior of Q(t), we have considered therefore 
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the "typical" average Q(t)t yp = exp((ln Q(t)}, where the brackets denote an 
average over many sets of N elements. Note that an arithmetic average over 
M sets of N elements can not be employed here, since it leads to a result 
identical for a larger system with M x N elements (see [|J). For a discussion 
of typical averages see 
the following. 



26 . For simplicity, we shall drop the index "typ" in 




Fig. 1. Plot of —\n[Q(t) /Q(0)] as a function of t in a double logarithmic presentation 
for (a) the trapping model in d = 1 and d — 2, and (b) the hierarchical constraint 
model, for several system sizes. For the trapping model, the system sizes are N = 
2 ■ 10 3 (open square), 2 ■ 10 5 (open circle), 2 • 10 7 (open up triangle), 2 • 10 9 (open 
down triangle) in d = 1, and N = 9 ■ 10 2 (full square), 9 • 10 4 (full circle), 9 ■ 10 B 
(full up triangle) in d — 2. For the hierarchical model, the system sizes are N — 10 2 
(full square), 10 3 (full circle), 10 4 (full up triangle), 10 5 (full down triangle). 



Figure 1 shows — ln[Q(t)/Q(0)] as a function of t in a double logarithmic 
plot for (i) the trapping model in d = 1 and d = 2, and (iii) the hierarchi- 
cal constraint model, both for several system sizes. In all cases, a crossover 
from an exponent (3 < 1 (at small t) towards (3 = 1 (at large t) can be eas- 
ily recognized. The crossover time t x shifts towards larger values when N 
increases. 

To study the crossover behavior in a more quantitative manner, we have 
plotted in Fig. |^ the local exponents (3 obtained from the local slopes of Fig. [l], 
as a function of t. In both systems, for a fixed system size N, (3 first decreases 
with t, reaches a minimum value at a certain time that can be identified with 
t x , and then increases monotonically with time towards (3 — 1. The figure 
shows that the minimum value of (3 has not yet reached its asymptotic value 
predicted for infinite systems, i.e., (3 — 1/3 (d = 1) and (3 = 1/2 (d = 2) for 
the trapping system and (3 = 1/2 for the hierarchical system. 
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Fig. 2. Plot of the local exponents (5 calculated from the successive slopes of the 
corresponding curves in ([j]) , (a) for the trapping model and (b) for the hierarchical 
model. The horizontal dashed lines represent the corresponding asymtotic (TV — » 
oo, t — > oo) values of j3. 



To show the dependence of the crossover time t x on the system size TV 
we have plotted, in Fig. [| the values of t"/( Q+7 ) as a function of In TV. The 
crossover time was obtained numerically from the position of the minima of 
the curves in Fig.H The resulting straight lines are in full agreement with 
the prediction of (0), supporting our analytical approach. 
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In the following we discuss the relevance of our results to Monte-Carlo 
simulations and experiments. There exists a long standing puzzle in Monte- 
Carlo simulations of the trapping problem in d = 2 and 3, that the predicted 
stretched exponential could not be observed || , even for survival probabilities 
Q(t)/Q(0) down to 1(T 21 in d = 2 [8b] and 1CT 67 in d = 3 [8g]. 

The finding of the logarithmic dependence of Q(t) on the system size 
N explains this puzzle. The Monte-Carlo simulations in d = 2 and 3 were 
typically performed on 10 3 configurations with about 10 4 traps, which is 
equivalent to having a single system with N ~ 10 7 traps. Using ([Toj) , we 
expect for N = 10 7 traps Q(t x )/Q(0) — 10~ 7 in d = 2. Indeed, for times 
above t x the exponent (3 aproaches unity as predicted by our theory and 
seen clearly in Fig. ||a. Moreover, for this system size j3 never reaches the 
predicted thermodynamic value f3 = 0.5, the minimum value of (3 is about 
0.65. For d = 3, Q(t x ) /Q(0) = 10~ n thus for smaller survival values (t > t x ) 
one again expects increasing values of (3 approaching unity. This explains the 
exponential decay found in the early Monte-Carlo simulations. Our results 
show that this is not an artefact but due to the finite size of the system. 
Moreover, they clearly indicate that the thermodynamic limit can not even 
be reached in one-dimensional macroscopic systems. 

It would be of interest to test the above prediction experimentally by 
preparing experimental realizations where size effects can be controlled. Equa- 
tions (8) and (10) suggest that the behavior around the crossover can be 
measured experimentally. For the trapping problem in linear systems, which 
has been studied experimentally |27],^8| , we expect for 10 8 sites and concen- 
trations of traps c between 10~ 4 and 10~ 2 , that Q(t x )/Q(0) ~ 10~ 2 10~ 3 , 
which is a survival range that can be detected experimentally. For the non 
identical particles (case (iii)), we expect for 10 s sites and concentration of 
traps c between 10 -4 and 1 that Q(t x )/Q(0) ~ 10 _1 -=- 10~ 2 which is a sur- 
vival range that can be well detected experimentally. The same arguments are 
valid for the target problem and therefore a similar crossover from stretched 
exponential to exponential decay is expected in relaxation experiments in low 
dimensional geometries |29|| , Mesoscopic systems such as quantum dots, are 
also promising candidates for experiments where the crossover can be rele- 
vant. Identifying the logarithmic size dependence in experiments may provide 
support to the theories claiming that the observed stretched exponential is 
due to competing exponential processes, represented by (^). 

This work was supported by the German Israeli Foundation (GIF). 
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